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Abstract 

After a brief discussion of the computational complexity of Clifford 
algebras, we present a new basis for even Clifford algebra Cl{2m) that 
simplifies greatly the actual calculations and, without resorting to the 
conventional matrix isomorphism formulation, obtains the same com- 
plexity. In the last part we apply these results to the Clifford algebra 
formulation of the NP-complete problem of the maximum clique of a 
graph introduced in [3j. 

1 Introduction 

Recently Clifford algebras have been applied to many "hard" problems, see 
e.g. [9j that show that several NP-complete problems require only a poly- 
nomial number of Clifford products to be solved, or [S] that transforms the 
solution of the maximum clique problem of a graph into the solution of an 
equation in Clifford algebra. 

So understanding the computational complexity of an actual calcula- 
tion in Clifford algebras is of paramount importance both for practical and 
theoretical reasons. In what follows we start applying well known results to 
show that the number of real multiplications actually needed for the explicit 
evaluation of a Clifford product has well defined bounds. 

Subsequently we introduce a basis for even Clifford algebras and show 
that in this basis the calculation of a Clifford product requires the same 
number of multiplications needed in the product of isomorphic matrices thus 



achieving this almost optimal result and hitting the upper bound exhibited 
previously. This base is made only of pure spinors and could be interesting 
in its own right. 

The final part follows the path opened in [3j where the maximum clique 
problem was formulated as a Cartan equation in Cl{2m); here we take advan- 
tage of presented results to achieve a more general and simpler formulation 
showing that between graphs and Clifford algebras deep relations exist. 

2 Complexity of Clifford Algebras 

Given a finite-dimensional, unital, associative algebra A over a field K, its 
multiplicative complexity C{A) is defined [S] as the essential number of 
multiplications needed to calculate the multiplication map of A, which is 
the bilinear map A x A ^ A and this definition can be made independent 
of coordinates. 

A well known result pQ (somewhat sharpened in [2]) states that for simple 
algebras 

C(^) > 2dimA- 1. (1) 

A Clifford algebra Cl{n) (see e.g. 0) is a simple algebra with n gen- 
erators 7i, 72, . . . , 7„, of dimension 2" and bound ([I]) implies that, given 
a,b £ Cl{n), the calculation of ab requires at least 2""*"^ — 1 multiplications, 
independently of the basis. We now compare two standard ways to calculate 
ab with this bound. 

The first possibility is to use a base {(} in which the generic Cl{n) 
element a is represented as: 

a = ^ aiCi where = 7i (2) 

where i € 21"] is a subset of [n\ = {1, 2, . . . , n} used as a multi-index and the 
2" "coordinates" ai € K. The actual computation of the product ab asks 
for the algebra multiplication table 

that can be easily computed from the generator properties 

IHj + lili ■= {lulj} = ■ 

This "direct" way to calculate the Clifford product ab is, on one hand 
intuitive and simple but, on the other, impracticable in all but the simplest 
cases given that all 2^" products QCj = Y[i& 7« T^jej 7i non-zero and all 
the 2^" multiplications of coordinates aibj have to be actually calculated. 
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From now on we consider Clifford algebras with even n := 2m, over field 
K and with vector space K'^"^. Even if the results that follow hold both for 
K = C and IR with signature 

72i-i = 1 72i = -1 i = l,...,m 

we will mainly address the real case, leaving to the reader the simple ad- 
justments for the complex case; given the IR^"* signature we will indicate 
the Clifford algebra with Cl{m,m). Since Cl{m,m) is graded isomorphic 
to K{2"^), the algebra of matrices of size 2™ x 2™ of dimension 2^"^, we 
can find again the lower bound to the multiplicative complexity of Cl{m, m) 
applying the bound (HD to K{2"'). 

The second possibility to calculate ab is to exploit this isomorphism and 
since standard matrix multiplication algorithms require O(r^) multiplica- 
tions for the calculation of the product of two r x r matrices it requires 
0(2^™) multiplications when applied to matrices that are isomorphic to 
Cl{m,mj^ a substantial complexity reduction with respect to the direct 
calculation, even if still a long way from the lower bound 

Anyway we can conclude that any actual calculation performed by means 
of a Clifford algebra formulation (see e.g. or [5) is sandwiched between 
these lower and upper bounds, respectively: 

2 22"^-! and 0(2^'"). 

3 Actual calculations in Cl{m,m): the Extended 
Fock Basis 

We have just seen that to perform actual calculations in Cl{m,m) the best 
is to take advantage from its (graded) isomorphism to matrix algebra K{2'^) 
but this is not always the case of choice fS] also because it's rather cumber- 
some. 

We propose here a method that allows to take the better of both worlds: 
on one side achieves the affordable result of 2^™ multiplications while, at 
the same time, maintaining the crisp formulation of 7 products. 

This can be achieved by means of a change of basis in Cl{m,m) that, 
exploiting the properties of Clifford algebras, produces pleasant properties 
as far as actual calculations are to be performed. This basis essentially 
extends to the entire algebra the Fock basis |4J of its spinorial part. 

""^we just mention here that there are faster matrix multiplication algorithms (see e.g. 
[10] and [7]) that, even if particularly well suited to matrices isomorphic to Clifford algebras 
that have sizes that are powers of 2, do not change substantially the general picture, in 
particular 0(2^"') = 0(8™) can be reduced to 0(7'") 
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We start defining the null, or Witt, basis of the vectorial part K"^^ of 
Cl{m,m) that takes the form: 



1 . 

with the properties 



and 



1 



(7: 



'2i-l 



72i) 



1,2,..., m (4) 



{qi,Qj} = and 



(5) 



that imply = qf = 0, at the origin of the name "null" given to this basis. 
With this basis i^^™ is easily seen to be the direct sum of two maximal 
Totally Null Planes (TNP) P and Q spanned by null vectors {pi} and {qi] 
respectively: 

since P H Q = {0} each vector v G K'^'^ may be expressed in the form 

m 

V = Yl {ociPi + PiQi) with ai,Pie K. 

i=l 

We now define the Extended Fock Basis (EFB) of Cl{m, m) to be given 
by all possible sequences 



'01'02 • ■ • V'n 



ipi e {qiPi,Piqi,Pi,qi} 



l,...,m 



and since every tpi can take just 4 values the basis contains 4™ = 2^"^ 
elements. Moreover we define as "signature" of an EFB element the vector 
(si, S2, . . . , Sjn) & {±1}"* where Si is the parity of ipi under the main algebra 
automorphism 7^ — 7j. 

We start with the simple example of CZ(1, 1) where the 4 EFB elements 
take the simple form {qiPi,Piqi,Pi,qi} and, with 7172 := 712 and standard 
matrix formalism, we can write 
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where the transformation matrix H between the standard 7 basis and EFB 

ij) can be written as I2 ® Hi where I2 is the identity matrix of size 2, H\ is 

11.111 

and has the properties 



the Hadamard matrix 



H = H^ 



14. 



We observe that the transformation H is orthogonal and thus EFB is 
a proper basis in Ci(l, 1) and its block structure allows to see the algebra 
as a direct sum of its even {1,712} and odd {71,72} parts that, in EFB, 
are mapped to the direct sum of the subspaces with +1 and —1 signature. 
For the general EFB in CZ(m,m) holds the following result, proved in the 
Appendix. 
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Table 1: Multiplication table of EFB elements in CZ(1, 1) 

Proposition 1 The EFB of Cl{m,m) is obtained from the standard basis 
'Pml by means of a matrix where: 

m 

H = l2m® Hm Hm = ®Hi (6) 

for which H = and ^HH = l22m. Vm is a permutation matrix de- 
fined recursively from Vm-i ('Pi = 14) and P23, the permutation matrix 
corresponding to permutation {1,3,2,4} and is given by: 

m— 1 

Vm = l2^ [Vm~l{ O P23)] ® I2 . 

The EFB is the direct sum of its 2™ subspaces with equal signatures each 
of them being the image of one of the 2™ matrices Hm appearing along the 
diagonal of the transformation matrix H . 

The following useful propositions are also simple to prove: 

Proposition 2 Given an EFB element ^ = ipi^2' ■ ■ ipm with signature 
(si, S2, • • • , Sm) and given another (not necessarily different) signature (ri, r2, . . . , r^) 
there exists one, and only one, EFB element <I> of signature (ri,r2, . . . ,rm) 
such that 7^ 0, moreover, in this case, the product is an EFB element of 
signature (siri, S2r2, . . . , Smrm)- 

In a generic Clifford product of EFB elements "if^ = ipiip2 ■ ■ ■ tpm 4>i4^2 - ■ ■ 4>m 
for i ^ j we have ipicpj = ztcpjipi and so 

• • • (t>l4>2 ■ ■■(f>m = ±V'l</'l'02</'2 • ' ' ^m'/'m 

and the only relevant products are thus ipicpi whose results appear in Table [H 
Prom the table it's easy to see that, given ipi of signature Sj there exist one, 
and only one, (pi of given signature rj such that 7^ and in this case 
the signature of the product is Sjrj. □ 

Given that Cl{m, m) is the direct sum of its 2™ subspaces with different 
EFB signatures one obtains: 

Corollary 3 Given an EFB element ^ = ipiip2 - ■ ■ ipm its Glifford products 
with any other EFB element <I) = 0x02 • • • (pm is not zero only for 2™ of the 
22™ elements ^ of the EFB. 
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Corollary 4 Given any EFB element ^ = "01 ^'2 ■ ■ ■ of Cl{m, m) either 
$2 _ vj/ (jj, xj,2 _ Q_ particular '^'^ = for all and only the 2™ elements 
of signature (1,1,..., 1) that thus form an even, Ahelian, sub-algebra of 
Cl{m,m). All other 2™(2™ - 1) EFB elements are null, ^'^ = q. 

This is a particular case of Proposition [2] for which we have shown that the 
only relevant products are ^l^i(f)i whose results appear in Table [TJ Here we 
are interested in ^j^, the diagonal of the table, and ipi'^ € {0, ?/'«}• Therefore 
since ^ = tp2 ■ ■ ■ V^m only if, for all i, tpi ^ we have proved 
that G {0, ±^}. To rule out the case ^'^ = we observe that if all 
ipi £ {QiPi^PiQi} then tpi^pj = ipjtpi for all i,j so that G {0, ^'}; the other 
parts of the proposition are trivial. □ 

An immediate consequence of these results is that when one wants to 
calculate the Clifford product ab in Cl{m,m), with o and b expressed in 
the EFB, for each of the 2^"^ coordinates of a we will need to calculate 2"^ 
multiplications that sum to just 2^*" for the entire Clifford product ab. 

In other words the multiplication table of the algebra (that is a table of 
size 2^"^ X 2^™, i.e. with 2^™ elements) has just 2^"* non zero elements. It 
shouldn't be difficult to prove that this is the minimum number of elements 
that are in general different from zero for any linear, invertible, transforma- 
tion of the standard 7 basis (wouldn't this be true one could get immediately, 
via matrix algebra isomorphism, an algorithm for matrix multiplication re- 
quiring less than the standard 0(2'^™') multiplications). 

So in Clifford products in EFB we can achieve the speed of matrix mul- 
tiplication without resorting to matrix isomorphism. 

Proposition 5 The 2^"* elements of EFB are simple (also: pure) spinors. 

We show first that all EFB elements \I' = V'iV'2 ■ ■ ■ V'm are Weyl spinors, i.e. 
defining F := 7172 • • -72™, that F^ = zt^*. We first note that 721-1722 = 
qiPi-PiQi and thus F = (qipi - piqi){q2P2 - P2q2) ■ ■ ■ (qmPm - Pmqm) ■ Then 
we note that {qiPi—piqi)ijjj = ipj{qiPi—Piqi) for i ^ j and consequently, that, 
as in previous proofs, only the products {qiPi — Piqiji^i are relevant. Since, 
depending on the values of -01! {qiPi — Piqi)tpi = iV'i this easily shows that 
F^ = ib^*. To prove now that the Weyl spinor ^ is simple it is sufficient 
to show that its associated TNP is maximal, i.e. of dimension m. For any 
"if = ipi'ip2 • • • 1pm let's call Vi the first null vector appearing in ^pi (thus Vi = qi 
for ipi G {qiPi,qi} and Vi = pi for V'i G {Piqi.Pi}) then Span{vi,V2, ■ ■ ■ ,Vm) 
is a TNP of maximal dimension m and for any v G Span{vi,V2, ■ ■ ■ , Vm) we 
have f\I' = 0, thus ^' is a simple spinor. □ 

We remark that, given the EFB element ^ and its associated maximal 
TNP Span(vi,V2, ■ ■ ■ , Vm) '■= M{^), there are in all 2"^ EFB elements whose 
TNP is M(^), that thus correspond to the same spinor ^ [Ij and that can 
be derived from ^ replacing every ipi with its counterpart with same first 
null vector Vi and opposite signature, i.e. pi piQi and qi qiPi. 
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We finally observe that it's immediate to change a base element of the 
standard 7 base, e.g. 7j, 7j, • • • , 7a: to a superposition of 2"* EFB elements 
substituting: 

• to each 72j-i the sum (pi + qi), 

• to each 72^ the sum (pi — qi), and 

• to each 7/ not appearing explicitly in 7^, 7^, . . . , 7^, the sum {qiPi + 
PiQi) = 1 

so that, for example 

72i-i72i = {qipi + piqi)iq2P2 + ^2^2) • • • (pi + - g'i) • • • {qmPm + Pmgm) 

and the product expands in a sum of precisely 2™ EFB elements all with 
the same signature. Viceversa every EFB element can be transformed in 
a linear superposition of 2"^ 7's by means of these properties clearly 
descend from the form of H 

We conclude observing that in Cl{m,m) the standard 7 basis and EFB 
have complementary properties. On one side in 7 basis the algebra is a direct 
sum of its m + 1 grades {K, vectors and multivectors) and all products of 
basis elements are non zero, on the contrary in EFB the algebra is a direct 
sum of 2™ parts with different signatures while the overwhelming majority 
of products of EFB elements is zero (only 1 into 2™" is non zero). 

4 A new formulation of the maximum clique prob- 
lem in Clifford algebra 

We start with a brief remind of the maximum clique problem of a graph and 
its formulation in Cl{m, m) appeared in [3] to which the reader is addressed 
for further details. Since a clique of a graph is a maximum independent set 
of its complementary graph, and this last representation is better suited for 
null vectors geometry, we will stick to it. 

Given a graph with m vertices and its adjacency matrix A with aij G 
{0, 1} one can define m vectors of Cl{m,m) 



and, given the properties ([H) of pi and qj G i^^™, one has: {zi,Zj} = Uij. 
With these vectors the maximum independent set of the graph A corresponds 
to the largest subset of vectors Zi that span a TNP in K'^^. 



m 




(7) 
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In the quoted paper we have shown that any maximal independent seel of 
A defines uniquely a maximal TNP plane in thus for example, if the set 
of vertices ji,j2, ■■■■,3k defines a maximal independent set then, indicating 
with Pi, . . . ,qj, . . . all Pi and qj that appear in any of the Zj-^ , zj^, ■ ■ ■ , Zj^ , 
Span{pi, . . . , Qj , . . .) is a maximal TNP i.e. of dimension m. Since in turn 
[1] the maximal TNP uniquely identify simple spinors of the spinor space S 
by means of the Cartan equation 

v^ = ve K^™, (pe S, v^cp^O (8) 

we have established an injective application Zi uiZi from maximal inde- 
pendent sets of a graph to simple spinors of Cl{m,m), i.e. elements of the 
Fock basis and thus of the EFB. Since a graph A is uniquely identified by 
the set of its maximal independent sets it follows that each graph uniquely 
determines a, usually not simple, spinor ^{A) € S that we may symbolically 
write as 

^{A) = ^ojz, 

I 

where the sum over / is extended to the set of maximal independent sets of 
A. We observe that this sum is not calculable in practice since it needs the 
set of all maximal independent sets, knowing which is equivalent to solving 
the maximum independent set problem. 

Let's consider the set of m bi-vectors piQi := Cj that we call Cj to alleviate 
the notation, it's easy to observe that they are the generators of an Abelian 
subalgebra of Cl{m,m) with the following properties: 

and 

CiQi = piCi = qiCi = Qi CiPi = Pi 

In correspondence to the Zi vectors ([7]) we define m multivectors o 

m 

where we assume = 1 and thus appears explicitly in Oj only if a^j = 1. 
We now prove the following: 

Proposition 6 o^Oj = for any i = 1, . . . ,m; for i ^ j OiOj = if, and 
only if ttij = 1. 

^an independent set is maximal if no further vertex can be added to it, the maximum 
independent set is the largest maximal independent set. 
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Table 2: Possible cases of products of multivectors Oj 



Since g^gj = and qi commutes with all other elements of o, it follows that 
OjOj = 0. If aij = 1 then ej appears in Oj and we can shift it to the right until 
it reaches qj and ejqj = 0. On the contrary let us suppose that OjOj / 0, 
then, given the properties of e bivectors, we deduce that ej doesn't appear 
in Oj which implies CLij = 0. □ 

This result follows from the asymmetry of the Clifford products: qiCi = qi 
whereas eiqi = and holds in general for any set of {oj} multivectors defined 
from any 0, 1 square matrix and not only for adjacency matrices of graphs. 
When the product OiOj ^ we use the notation OiOj = qiqje ■ ■ ■ e where all 
e's are shifted to the right and they all appear with power 1 (remember 
are idempotent); neither nor ej may appear in the set e • • • e if the product 
is not zero; Table [2] resumes the 4 possible cases that may occur. 

We introduce now a slightly different matrix A' that is essentially the 
starting adjacency matrix A with the lower triangle elements all set to 1, 
more precisely 

a'ij = aij for j > i and a^j = 1 for i > j 

and with this matrix we redefine the m multivectors that thus take the 
form Oi = 6162 • • • Ci-i qi e""J"^ • • • e^™-, in this fashion every graph A defines 
uniquely an element of Cl{m,m) given by 

m 

= ^0. 

i=l 

and we observe that each Oj can be easily written in the EFB, it suffices to 
substitute to each aij = 

e^y =1 = {qipi +piqi) 

and thus, if Oj has / indexes ji, j2) ■ ■ ■ ,ji such that aij = then Oj will be 
written as a sum of 2' EFB elements. We now use O to prove 

Proposition 7 The graph A has an independent set of order k<mif, and 
only if, 

m 
i=l 
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Let the independent set be identified by vertices ji, j2) ■ ■ ■ ,jk with ji < j2 < 
• • • < ife then contains at least the term Oj^ oj^ ■ ■ ■ Oj^ . // the independent 
set is maximal Oj^ Oj^ ■ ■ ■ oj^ is a simple spinor. 

We proceed by induction: let's start by = 2, in this case 

m m 
= C^Oif = "^of + ^{Oi,Oj} = OiOj 
i=l i=l j>i j>i 

and the sum contains only terms that refer to independent sets of size 2, i.e. 
links. Let's now suppose that the relation is true for an independent set of 
size k, i.e. 

= ■■■ + Oj,Oj^---Oj^ + ■■■ 

with indexes Ji < J2 < ■ ■ ■ < jk that identify an independent set of size k. 
We can then write 



i=l 



and let's consider the generic resulting term Oj^Oj^ ■ ■ ■ Oj^Oi, by previous 
corollary we know that Oj^Oi ^ only if Oj,,j = that, in our case, im- 
plies i > jk- This term can be written 



qj^e---e qj^e-'-e ••• g-,-^e---e qie---e 

and eji_ doesn't appear in qj^^e - ■ ■ e since ajj., = but the overall product is 
not zero only if ej^, is missing also from all the preceding o's. This means 
that also ji,j2, ■ ■ ■ jjk^i form an independent set with ji < j2 < • • • < jk < i 
that concludes the induction argument. The viceversa is trivial and this 
proves the first part of the proposition. 

To prove the second part we observe that if the independent set ji < j2 < 
■ ■ ■ < jfc is maximal no other vertex can be added to it, or, equivalently, no Oj 
can be inserted (in its proper position) to the product Oj^Oj^ ■ ■ ■ Oj^. without 
sending it to zero. Since 

OhOj.,---Oj^=qj^e---e qj^e-'-e ■■■ qj^e---e 

where in the set of the indices of g's and e's there all elements 1,2,... , m 
whether among the g's or among the e's. On the contrary if, e.g., i were miss- 
ing this would imply i > jk-, since Oj^, contains surely at least eie2 • • • e.j^-i 
and we could thus append Oj to Oj^Oj^ ■ ■ ■ Oj^ without zeroing it against the 
hypothesis of a maximal independent set. We conclude that Oj^Oj^ ■ ■ ■ Oj^ is 
an element of EFB that is a simple spinor by Proposition [5j □ 

We remark that the second part of the proposition gives a necessary 
condition that is not sufficient since there are non maximal independent 
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sets, e.g. Om, that are EFB elements. Another interesting observation is 
that if = necessarily all the elements of the expansion of are zero 
since, being all these terms linearly independent (because of the qi), no 
cancellation can occur among them. A more general consequence is: 

Corollary 8 Every nonzero simple spinor (p of Cl{m,m) can be written 
(j) = for some k. 

Let (p = Q1Q2 • ■ ■ Qk • • • em be an EFB element, then all of those O 

generated by graphs that have 1, 2, . . . , fc as their unique maximum indepen- 
dent set of size k, satisfy O'^ = cj). n 

This last observation shows that the relation between graphs and simple 
spinors can be much deeper than previously thought. 

A Appendix: proof of Proposition [1] 

To prove Proposition [T] we build constructively the general EFB in C/(m, m) 
using the isomorphism 

Cl{m + 1, m + 1) ^ Cl{m, m) (g) CI{1, 1) 

and so the standard 7 basis in C/(2,2) is given by 7 (8) 7 and using the 
relations found for C/(l, 1) 

V'OV' = ^{H'j)(E){Hj) = ^{H(g)H){-f(g)j) = ^(l2®i/i®l2«'^^i)(7«'7) (9) 

and it's easy to see, calling P23 the symmetric permutation matrix corre- 
sponding to permutation {1, 3, 2, 4}, that P23(l2 ® ^^i)-P23 = -ffi <8)l2 so that 
we may write for the new transformation matrix H ® H = [-^23(12 ®> 
-^^1)^23] ® Hi and, with easy passages, 

H(^H = {t2®P2^®t2){t2m2®Hi®Hi){l2®P2^®t2) := V2{U®H2)V2 

where V2 is a symmetric permutation matrix and H2 = Hi ® Hi is the 
Hadamard matrix of size 2^. We can left multiply ([9]) by P2 to get 

V^{il, ® ^/;) = i(l4 H2) ^2(7 ® 7) 

where 

V2{l®l)=V2 (1,712,71)72)'^ <^ (1,734,73,74)'^ = 

(1,734,712,71234, 73,74,7123,7124, 7l , 7l34 , 72 , 7234 , 713,714,723,724)'^ • 
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It is simple to use this recursive construction to build any base in Cl{m, m) 
thus proving the first part of the proposition. To calculate explicitly Vm we 
first note that, assuming Vi = I4, we may write in Cl{m, m) 

rn 1 m m 

and 

m m, ^ 

<^H = (8)(l2 <8) = I2 <8) [ (8) {Hi (g) I2)] <8) Hi 

and that 

m—1 m—1 m—1 m—1 m—1 

(g) {Hi ® I2) = ® (P23I2 (g i^l-P23) = ( ^ P23)[ ® (I2 ® Hi)]{ «) P23) 
m—1 rr^ m—1 

= ( ® P23)7'^_i(l2— 1 ^^m-OT'm-ll ® ^^23) 

so that 

m m—1 

(8)-ff = I2 (8) [ ® (i^i (g I2)] » i^i 

= I2 [( (8 P23)7'^-l(l2— 1 <8 H^_i)Vm-l{ ® P23)\ ® Hi 

= {I2 ®[{® P2z)Vl,_i] ® I2KI2- ® H^){l2 [Vm-l{ ® -P23)] O I2} 

= {[I2 ® C®P2Z) ® l2][l2 ® Vm-l (8) 12]}(12- O ^f'm) 

m—1 

{[I2 0:Pm-l ® l2][l2 8) ( O ^^23) O I2]} 
:= P^(l2- (8 i/m)7'm 

that provides the desired definition of Vm. To prove the last part of the 

m 

proposition one just need to show that the 2™ subspaces obtained by Vm{®i^) 
all have the same signature and one can proceed by induction. The proposi- 
tion has been shown true for CI{1, 1) so if it's true for Cl{m, rn) this means 

m m+1 

that Vm{®'tp) satisfy the condition. To show that it holds for ® tp one can 

m+l m m 

use the associativity of the external product ® = ® {®'^) = ((8)'0) "0 

m+l 

to argue that in Vm+i{ ® ^) both the first m elements and the last m must 
have the same signature; it follows that all m + 1 elements have the same 
signature and the property holds. □ 
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